DMRG Study of Critical Behavior of the Spin-1/2 Alternating Heisenberg Chain 
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We investigate the critical behavior of the 5 = 1/2 alternating Heisenberg chain using the density 
matrix renormalization group (DMRG). The ground-state energy per spin eo and singlet-triplet 
energy gap A are determined for a range of alternations 8. Our results for the approach of eo to 
the uniform chain limit are well described by a power law with exponent p w 1.45. The singlet- 
triplet gap is also well described by a power law, with a critical exponent of p w 0.73, half of the 
eo exponent. The renormalization group predictions of power laws with logarithmic corrections can 
also accurately describe our data provided that a surprisingly large scale parameter So is present in 
the logarithm. 

PACS numbers: 75.10.Jm, 75.40.-s, 75.40.Cx 
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I. INTRODUCTION 

The approach of quantum spin systems to criticality is 
an interesting and rather underexplored problem in com- 
putational physics. In this paper we consider the critical 
behavior of one of the simplest quasi- ID quantum spin 
systems, the alternating Heisenberg chain with spin-1/2, 
as it approaches the uniform chain limit. This model is 
of special interest for the study of critical behavior, since 
much is known analytically about the gapless uniform- 
chain limit through the Bethe Ansatz. This system is also 
important for studies of the magnetic spin-Peierls effect, 
and gives a reasonably accurate description of magnetic 
interactions in many dimerized quasi-lD antiferromag- 
nets. 

The alternating Heisenberg antiferromagnet is a simple 
generalization of the uniform Heisenberg chain, with the 
nearest-neighbor exchange constant J replaced by two 
alternating values. The spin-1/2 system is defined by 
the Hamiltonian 

L/2 

H = - J Y, ((! + 5 ) ' & + (!-*) & ■ 3x+i 

i=l 

(1) 

and has a spin-singlet ground state with energy Eq and a 
nonzero singlet-triplet gap E\ — Eq for any alternation!]] 
< 8 < 1. 

For our study we introduce a finite-lattice, scaled 
ground-state energy per spin, 



(2) 



(3) 



The tilde indicates that these energies are scaled by J 
rather than J = (1 + 5) J, which would divide our results 



e (L,5) = E (L,5) / LJ 
and a corresponding singlet-triplet gap 

ACM) = (E 1 (L,5)-E (L,5)) I J 



by a factor of (1 + ft) (albeit giving the same leading crit- 
ical behavior). In addition, energies quoted with a single 
real argument and a tilde (e.g. &o(S)) are bulk limits of 
the scaled alternating chain, those with no tilde and a 
single integer argument (e.g. eo(L)) are finite-!/ uniform 
chain results, and those with no arguments (e.g. eo) are 
bulk limits of the uniform chain. 

The uniform spin-1/2 Heisenberg antiferromagnetic 
chain, which we recover at S = 0, is the best under- 
stood ID critical quantum spin system. It has a ground- 
state energy per spin of eo = 1/4 — In 2 and a band of 
gapless spin-triplet excitations with dispersion relation 
uj(k)/ J = 7r/2 I sin(fc)|. The approach of noncritical mod- 
els to this limiting case is not well established, and has 
been the subject of surprisingly few theoretical and nu- 
merical studies. 

Early analytical studies of the effect of a small alterna- 
tion Sjfn the uniform chain were reported by Cross and 
Fishcra and Black and Emerya. Cross and Fisher used a 
Jordan- Wigner transformation to map the original spin 
problem onto a pseudo-fermion Hamiltonian, and ap- 
proximated the latter by the exactly solvable Luttinger- 
Tomanaga model. This approach, which unfortunately 
involves uncontrolled approximations, yields critical ex- 
ponents (defined by f(S) oc 8 P ) of p = 4/3 for the ground- 
state energy and p 2/3 for the singlet-triplet gap. 

Black and Emerya related the critical behavior of the 
alternating Heisenberg chain to the 4-state Potts model, 
and found logarithmic corrections to these power laws, 



e -e (5)^S 4 / 3 /\logS\ 
for the ground-state energy per spin, and 
A(S) oc(5 2 / 3 /|k>g<$| 1/2 



(4) 



(5) 



for the gap. Note that at this order the gap scales as 
the square root of the ground-state energy defect. More 
recent theoretical work by Affleck et ala has shown that 



2 



Eqs. are leading-order predictions of the renormal- 
ization group. The overall constants in these results and 
implicit in the logarithms are non-universal, and have not 
yet been determined analytically for this model. 

Several previous numerical studies have investigated 
the critical behavior of spin-1/2 alternating Heisenberg 
chains using bulk-limit extrapolations of exact diagonal- 
ization results on. systems up to about L — 30 in extent 
(see Barnes et ala and Yu and Haaaj and references cited 
therein). We shall see that important systematic errors 
can arise from extrapolations using these relatively small 
systems, for example in estimates of critical exponents. 

Studies of much larger spin-1/2 alternating Heisenberg 
chains have also been published using the DMRG algo- 
rithm, although there has been little systematic study 
of the critical behavior of the simple alternating chain 
model of Eq.(l) using DMRG. The single published ref- 
erence on this topic is the work of Uhrig et ala, who 
estimate a gap critical exponent of 0.65. Their numerical 
energies however deviate systematically from this power 
law at small alternation 8 (see Fig. 3 of Uhrig et ala). 
In related work, Chitra et alB used the DMRG method 
to study the effects of dimerization and frustration on 
a generalized alternating chain model with next-nearest- 
neighbor couplings, and Lou et alx3 studied the gap in- 
duced by a staggered magnetic field. 

In this work we present a systematic DMRG study of 
the critical behavior of the original alternating spin-1/2 
Heisenberg chain of Eq.(l), and compare the theoretical 
predictions Eqs. (4,5) to numerical results for the ground- 
state energy per spin and singlet-triplet gap of this model, 
using systems up to L = 192 in extent. 



II. NUMERICAL METHOD 

For our numerical study we employed the density ma- 
trix renormalization group (DMRG) algorithmd, which 
is a very effective method for studying critical behavior 
in quasi-lD quantum spin systems. Numerical determi- 
nation of the critical behavior of the alternating chain, 
or any similar quantum spin system, is a daunting com- 
putational task. One must accurately determine energy 
eigenvalues on quite large systems, since characteristic 
lengths typically diverge at critical points. An extrapo- 
lation through a series of fixed- Z/ results is then required 
on sufficiently large lattices to insure that one is in an 
asymptotic regime in which finite size effects can be accu- 
rately parametrized and eliminated. The computer mem- 
ory requirements for diagonalizing these large systems are 
such that the detailed critical behavior of relatively few 
quantum spin systems has been explored numerically. 

The DMRG algorithm has previously been applied 
to various interacting fermion systems, including one- 
dimensiorpJ spin chainsu j— lattice modelsEj, quasilinear 
moleculesEJ and nucleicMj. The essential concept in 
DMRG is to "grow" a small, finite system into a larger 
one by the iterative incorporation of new lattice sites. At 



each such iteration one retains only the m most relevant 
basis states for spanning the targeted energy eigenstate. 
(These basis states are chosen according to a density- 
matrix weight.) This selective sampling of Hilbert space 
yields accurate energy eigenvalues on systems which are 
well beyond the limits of exact diagonalization. For de- 
tails of the DMRG method we refer the reader to the 
original papers of WhiteEl and to-a series of lectures re- 
cently compiled by Peschel et alLJ 

Our numerical implementation proceeds as follows. We 
divide the spin chain into blocks A — a — B — b, where A 
and B denote the "system" and "environment" blocks, 
and a and b are elementary blocks to be added to A and 
B respectively. Blocks b and A are linked by periodic 
boundary conditions. We take spin dimers (two lattice 
sites of spin s = 1/2 each) as our elementary blocks, and 
use the infinite algorithm to grow the spin chain, while 
targeting the lowest-lying spin-0 or spin-1 state. In the 
case of the spin-0 ground state of the alternating Heisen- 
berg chain, we found that subsequent sweeps with the 
finite algorithm do not lead to much improvement in the 
results of the infinite algorithm. In contrast, for the spin- 
1 state we found that sweeps with the finite algorithm 
were important for convergence. Our DMRG implemen- 
tation uses largej-sparse matrices, so the sparse matrix 
package ArpackEj was employed. The time consuming 
matrix- vector multiplications and the Arnoldi/Lanczos 
algorithm Arpack were parallelized. 

The ground-state energy per spin eo(L,8) of a spin 
chain of length L and alternation 8 is of special in- 
terest for our study of critical behavior (see Eq. @). 
DMRG yields approximate values eo(m; L,8) that typi- 
cally converge exponentially fast from above as one in- 
creases the number m of states retainecoLIl. We also 
observed this behavior in the present study. We com- 
puted eo(m; L, 8) numerically using DMRG on chains of 
length L = 28, 48, 96, 144 and 192, and alternation 
8 = 2~ 3 , 2~ 4 , 2~ 5 , 2- 10 . At each L and 8 we increased 
m in steps of 10 starting at m = 30 until a fit of the form 
eo(m;L,8) = eo(L,8) + cicxp(— C2m) gave sufficiently 
stable coefficients for our desired accuracy; this m ex- 
trapolation yielded our DMRG energy estimate £q(L,8). 
We found that the maximum m needed for convergence to 
a given accuracy increases with increasing chain length 
L and decreasing alternation 8. For our extreme case 
L = 192 and 8 = 2~ 10 , adequate convergence was not 
achieved until m = 150, and we retained the maximum 
of m = 170 states in this case. This resulted in a sparse 
matrix problem of dimension w 9 x 10 4 at each DMRG 
iteration. Finally, we also confirmed recovery of exact 
L = 28 alternating chain results on allowing our DMRG 
code to iterate to the full Hilbert space. 



III. RESULTS 

Table | gives our DMRG results for the ground-state 
energy per spin for different alternations 8 and chain 
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FIG. 1: Approach of eo(I/,5) to the uniform-chain limit eo = 
1/4 — ln(2) with decreasing alternation S and increasing chain 
length L. The dashed line shows a fit to a power law, coS p , 
which gives Co = 0.39 and an exponent of 1.45. The full line 
shows the renormalization-group exponent 4/3 (displaced for 
presentation, without logarithmic corrections). 
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FIG. 2: Critical behavior of the singlet-triplet gap A(L,S) 
as a function of alternation S for spin chains of length L. 
The dashed line shows a power law fit, co8 p , which gives 
co = 1.94 and exponent p — 0.73. The full line shows the 
renormalization-group exponent p = 2/3 (displayed as in 
Fig. |). 



lengths L. For the rather large alternations 6 = 2~ 3 
and 2 -4 we confirmed convergence to accurately known 
energies with increasing basis state number m, as well as 
convergence with increasing chain length L. Finite size 
effects were more pronounced at smaller 6, as expected 
since the system is closer to criticality. We estimate that 
the energy errors in Table | are a few units in the last 
digit, based on the difference between the maximum-m 
DMRG result and the exponential m — > oo extrapolation 
which we quote. 

Fig. shows the approach of eo(L,S) to the limit 
eo = 1 /4 — In 2 as a function of alternation S for spin 
chains of increasing length, L = 28, 48, 96, 144 and 192. 
The asymptotic large- L envelope evident in this figure 
for S > 2~ 8 clearly shows the bulk limit. The dashed 
line is a power-law fit cqS p to this envelope, which gives 
Co = 0.39 and exponent p = 1.45. The renormalization- 
group prediction p = 4/3 from Eq. (|4|) (without the log- 
arithmic factor) is also shown, displaced for clarity of 
presentation. Clearly this gives an inferior description of 
the DMRG data over the range considered here. 

Fig. [2] shows the singlet-triplet gap A(L,S) as a func- 
tion of alternation S for spin chains of length L = 28, 48, 
96, 144 and 192. The envelope of these curves is the bulk 
limit A((5), which is evident for S > 2~ 7 . A power-law fit 
to A (6) as for eo(S) gives Co = 1.94 and p — 0.73, shown 
as a dashed line in Fig. ^. The renormalization-group pre- 
diction p = 2/3 from Eq.(||) is also shown, displaced for 
presentation. Evidently the renormalization-group expo- 
nent (without the logarithmic term) again gives a less 
accurate description of our DMRG data. Note that the 
two exponents obtained in our fits are related by a fac- 
tor of two, which implies that the eo(S) defect scales as 



A((5) 2 . This relation also follows from the leading-order 
renormahzation group, Eqs. (|JJ^). 

It is especially interesting to determine whether there 
is numerical evidence for logarithmic corrections to pure 
power-law behavior, as predicted by the renormahza- 
tion group in Eqs.(||J^). We first consider the bulk-limit 
ground-state energy per spin, predicted to asymptotically 
approach the uniform chain limit as 



e - e (<5) 



£4/3 



(0) 



Here we have introduced an overall constant a and scale 
parameter £0, which we will estimate from our DMRG 
data. On rearranging Eq. (j(^) we obtain the easily visu- 
alized form 



£4/3 



eo - e (S) 



In 5 — In £ 



(7) 



Fig. I shows our DMRG results for 6^ 3 /{e - e (5)) 
(the l.h.s. of Eq. (@)) versus In £. The data clearly dis- 
agree with the leading-order renormahzation group pre- 
diction Eq. @ over this range of £, since the points do not 
lie on a straight line. Assuming that the two smallest-£ 
points are close to asymptotic, we estimate a rs —2.2 and 
£0 « 110. A similar fit to the singlet-triplet gap gives 



£2/3 



U(£) 



( 



(8) 



with ctg a p Rj —19.4 and £ « 115. With these constants 
the theoretical renormalization-group results Eqs. (|^J^) 
are barely distinguishable from our power-law fits in 
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FIG. 3: Estimation of the constants a and So in Eq. 
our DMRG data. 



from 



Figs. and ||. There are no predictions of these con- 
stants in the literature to our knowledge. 

Eq. (H) can be reexpressed as an effective power 
Peff (5) = 4/3 + l/\n(S /S). With S = 110, the range 
S = 10~ 3 -> 1CT 1 corresponds to p eff = 1.42 -> 1.47, 
which may explain our good numerical agreement with a 
pure power of exponent p = 1.45. 

We have also investigated finite size effects in the uriL. 
form Heisenberg .chain (S — 0). Woynarovich and EckletEl 
and Affleck et aln quote Bethe-Ansatz predictions for the 
leading finite size contributions to the ground-state en- 
ergy per spin and singlet-triplet gap, 



eo 



and 



e (L) = 



12L 2 



l + (9((lnL)- 



(9) 



ex(L) - e Q = ^ (l - - (lni)- 1 + 0((ln L)- 2 )) . (10) 

Fig. || shows our DMRG data for these energy defects on 
uniform chains of length L = 28, 48, 96 and 144, together 
with the Bethe-Ansatz results. 

The agreement between the leading-order Bethe- 
Ansatz predictions for finite size corrections to the 
ground-state energy per spin (Eq. (||)) and our DMRG 
results is evidently very good on the shorter chains. The 
discrepancy evident at L = 144 may be due to conver- 
gence problems encountered by the DMRG algorithm 
when applied to this asymptotically gapless system on 
large lattices. The DMRG results for the spin-1 energy 
defect depart significantly from the leading-order term in 
Eq. ( |To| ) , but are consistent with this prediction when we 
include the 0(L _2 (ln L)~ l ) correction. 

Finally we consider the dependence of the ground- 
state energy per spin e (L,5) on the chain length L for 
fixed, nonzero alternation 5. This is especially relevant 
to exact diagonalization studies, which extrapolate to the 
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FIG. 4: Finite-size energy defects \e 3 (L) — eo| of the low- 
est spin-s levels on the uniform chain versus chain length L. 
The data points are DMRG (+:spin-0, x :spin-l) and the lines 
are Bethe-Ansatz results (leading order (LO}, solid; next-to- 
leading order (NLO), dashed), from Eqs.(§0). 



bulk limit from rather small systems of at most about 
30 spins. To test the accuracy of finite size extrapola- 
tions we used DMRG to compute the ground-state en- 
ergy per spin and singlet-triplet gap on chains of length 
L = 28, 32, 36, 40, 44 and 48, and fitted the results to 
f(L) = a + bexp (—L/c)/L and g(L) = a + fecxp (— L/c) 
respectively. (These forms were used to extrapolate ex- 
act diagonalization results r-to the bulk limit by Barnes 
et alB and Yu and Haas.H) Both functions yield rea- 
sonably good fits to our six data points in the range 
L = 28,..., 48. However, on comparing this extrapo- 
lation with our DMRG results on L = 96, 144 and 192 
lattices, we noted clear discrepancies as the alternation 
5 decreases. This is likely due to a rapid increase in the 
characteristic length (modeled by c in the exponents of 
f(L) and g(L)), which makes subleading terms in the 
asymptotic behavior more important. Our results thus 
suggest caution in attempting to establish critical behav- 
ior from studies of relatively small systems. 



IV. SUMMARY 

We have employed the DMRG algorithm to determine 
the ground-state energy per spin and singlet-triplet gap 
of the alternating Heisenberg chain, and studied the crit- 
ical behavior of this model in the limit of small alterna- 
tion. We find that the approach of the bulk-limit ground- 
state energy per spin and singlet-triplet gap to the uni- 
form chain limits are well described by power laws in the 
alternation parameter S over the range 0.008 < 6 0.1, 
and scale approximately as S 1A5 and 5° 73 respectively. 

The renormalization group predictions of power laws 
times logarithmic corrections also appear consistent with 
our results, provided that a surprisingly large scale factor 



5 



Sq is present in the logarithms. 
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e (L,5) 


ei(L,S) 
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28 


-0.4110961 


-0.3976678 


2-' A 


48 


-0.4110928 


-0.4032882 


2- 3 


96 


-0.4110928 


-0.407191 


2- 3 


144 


-0.4110928 


-0.4084913 


2- 3 


192 


-0.41109284 


-0.4091416 


2- 4 


28 


-0.4239186 


-0.4149855 


2- 4 


48 


-0.4238627 


-0.418895 


2- 4 


96 


-0.4238617 


-0.4213883 


2- 4 


144 


-0.4238619 


-0.422211 


2- 4 


192 


-0.4238618 


-0.422623 


2"° 


28 


-0.4325593 


-0.4258695 


2- 5 


48 


-0.4323071 


-0.429086 


2~ 5 


96 


-0.4322899 


-0.430754 


2- 5 


144 


-0.432290 


-0.431263 


2- 5 


192 


-0.4322906 


-0.431518 


2"° 


28 


-0.4378673 


-0.4319820 


2- 6 


48 


-0.437370 


-0.435005 


2- 6 


96 


-0.437284 


-0.436321 


2- 6 


144 


-0.437283 


-0.436652 


2- 6 


192 


-0.437285 


-0.43681 


2"' 


28 


-0.4408870 


-0.435219 


2- 7 


48 


-0.440254 


-0.438202 


2- 7 


96 


-0.440074 


-0.439410 


2~ 7 


144 


-0.440058 


-0.439664 


2~ 7 


192 


-0.440064 


-0.439768 


2-° 


28 


-0.4425056 


-0.4368847 


2- 8 


48 


-0.441826 


-0.439863 


2- 8 


96 


-0.44158 


-0.4410442 


2- 8 


144 


-0.441550 


-0.441275 


2- 8 


192 


-0.441550 


-0.441360 


2 -» 


28 


-0.4433439 


-0.4377293 


2 -9 


48 


-0.442650 


-0.440709 


2- 9 


96 


-0.442384 


-0.441885 


2 -9 


144 


-0.442340 


-0.442109 


2- 9 


192 


-0.442336 


-0.442190 


2 -iu 


28 


-0.4437703 


-0.4381546 


2 -io 


48 


-0.443073 


-0.441136 


2 -io 


96 


-0.44281 


-0.442311 


2 -io 


144 


-0.442756 


-0.442532 


2 -io 


192 


-0.442741 


-0.442614 





28 


-0.444201 


-0.4385820 





48 


-0.443504 


-0.441566 





96 


-0.44323 


-0.442741 





144 


-0.44318 


-0.442965 



TABLE I: 

Extrapolated DMRG results for the lowest energy per spin 
in the spin-s sector, e s (L, 8) = e s (L, S)/(l + 8), for various 
alternations 8 and chain lengths L. The estimated error is a 
few units in the final digit. 



